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Abstract 

Bennett's pebble game [1, 2] was introduced to obtain better time/space 
tradeoffs in the simulation of standard Turing machines by reversible ones. 
So far only upper bounds for the tradeoff based on the pebble game have 
been published. Here we give a recursion for the time optimal solution of 
the pebble game given a space bound. We analyze the recursion to obtain 
an explicit asymptotic expression for the best time-space product. 

1 Introduction 

In Bennett's pebble game, we have a board of n squares labeled si, . . . , s„. Each 
square may be empty or contain one pebble. A move in the game consists of 
adding or removing a pebble. A pebble can be added to or removed from square 
Sj if i = 1 or if there is a pebble on square Initially, the board is empty. 

The desired final configuration has exactly one pebble on the board on square 
Sn- For the motivation and application to reversible computing, we refer the 
reader to [1, 2]. 

The pebble game is trivial without resource constraints. At the i'th 
step for 1 < i < ra we can place a pebble on square Si. At the ra + i'th step for 
l<i<ra— Iwe can remove the pebble on square Sn-i- This solves the game 
in 2ra — 1 steps using n pebbles. The time resource T is the number of steps to 
solve the game, while the space resource S is the number of pebbles used, that 
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is the maximum number of pebbles on the board at any time. For any solution 
S of the pebble game, let T(S) be the number of steps and S(S) the number of 
pebbles used. Our goal is to determine 

F{n,S) = mm{T{S) \ S{S) < S). 

Levine and Sherman [2] expressed Bennett's [1] upper bound on F(n,S) by 
a recursion which they were able to solve. We will obtain an exact recursive 
expression for F(n, S) which can be evaluated explicitly for any desired n and S . 
An asymptotic expression for the minimum time space product is obtained. A 
consequence of this analysis is that an upper bound for the minimum time-space 
product of quantum factoring (see [3]) is 2^^(^"'"°(^))ra"^, where n is the number 
of bits in the number to be factored. We defer discussion of this observation to 
another note. 

2 Analysis of the pebble game 
2.1 Solutions of the pebble game 

To obtain a recursive expression for F(n, S) it is useful to model solutions S by 
representing the residence of each pebble on a square as an interval. Formally, 
we view 5 as a — 1 matrix, also denoted by S, of dimensions n X [— oo, +00]. 
Si J = 1 means that there is a pebble on square Si after step j. For 5 to be a 
valid solution, we need 

(i) Initial configuration: There exists Ff, such that for t < Ff,, Si^t = 0. 

(ii) Final configuration: There exists Te, such that for t > Fg, Si^t = ^n,i^ with 

the usual delta function. 

(iii) One move at a time: If Si^j 7^ ^i^j+i then for all i' ^ i, Si'j = Si'j^i. 

(iv) Add constraint: If Sij = and Sij-^i = 1, then either i = 1 or Si-ij = 1. 

(v) Remove constraint: If Sij = and Sij-i = 1, then either i = 1 or 
•Si-l,]-! = 1- 

The space requirement of S is given by the maximum number of I's in any 
column of S, while the time requirement is the number of inequalities between 
adjacent elements in the rows. 
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Consider square Si. The i'th row of S can be viewed as a collection of 
intervals of I's, representing the residence of a pebble on Si. Thus a (residence) 
interval [k,l] of Si satisfies (1) k < I, (2) Sij = 1 for k < j < I, (3) Si^k-i = 
and (4) 5,- ;+i = 0. 

Let J, j' > 0. If / = [A;, /] is an interval of Si and / = [k + j,l — j'] is an 
interval of Sj_i,then / does not contribute anything to the solution. A better 
solution in terms of time and space is obtained by removing / (by setting the 
corresponding entries of S to 0). This allows us to add another constraint: 

(vi) If / is an interval of Si and / is an interval of then J <^ I. 

In depicting solutions it is convenient to represent the intervals as horizontal 
lines at a vertical offset, where no two endpoints have the same horizontal offset. 
See Figure 2.1 for two examples. 

2.2 A recursion for F{n,S) 
Theorem 2.1 F(n,S) is given by 

for n = 1 and S > 1, 
for n > 2 and S = 1, 
m, S — 1) \ 1 < m < n) 
for n > 2 and S > 2. 

Proof. For ra = 1, the solution is to place a pebble on si in one step. For 
ra > 1, at least 2 pebbles are required. 

The recursive expression for F(n, S) is obtained by decomposing the 
solution matrix. The basic idea is shown in Figure 2.1. We begin by proving 
that F(n, S) > mm(F(m, S) + F(m, S - 1) + F(n - m, S - 1) \ 1 < m< n). Let 
n > 1 and S > 1. Let 5 be a time optimal solution of the pebble game with 
S(S) < S. Let tf, be the least t such that Si^t = 1- We can assume that for 
t > tf,, the fth column of S always has at least one 1. If not we can obtain a 
more time efficient solution by setting all the entries preceding that column to 
0. Let te be the largest t such that Si^t = 1- Due to the move constraints, for 
t > te, Si^t = ^n,i- Define m(S,t) = min(i | Si^t = 1) and let t^ be the maximum 
t with tfj < t < te for which m(S,t) achieves its maximum in this range. Let 
m = m(S,tm)- Observe that m < n, else we can obtain a better solution by 
setting Si^t = for i < m and t > t.^. We can extract three solutions of smaller 
problems from S. 



F(n,S) = I 
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min(F(m, S) + F(m, S - 1) + F(n 
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n=4, S=3: 



o 



O 



n=8, S=4: 



Let S[ be the submatrix of S consisting of the first m rows between 
columns 4 and ■ Let Si be the extension of S[ to a solution of the pebble 
game obtained by letting Si-i^t = for t < tf, and Si-i^t = ^m,i for t > t^. Then 
Si is a solution of the pebble game with m squares using at most S pebbles. 

Let S2 be the the submatrix oiS consisting of the first m rows between 
columns t.^ and tg. Let S2 be S2 with the columns in reverse order and extended 
to a solution as above for S[. Then S2 is a solution of the pebble game with m 
squares. We claim that it uses at most S — 1 pebbles. To see this we show that 
for t > tm, the maximum i(t) for which Si^t = 1 satisfies i > m. Suppose that 
this is not the case. Let ti and ^2 be the right endpoint and the left endpoint of 
the two intervals of si surrounding t^- Because of the move constraints, Si^t = 
for ti < t < t2 and i < m. This implies that the columns between ti and ^2 
have only one 1. Since S > 2, we can set Si^t to 1 for ti < t < t2. This reduces 
the number of steps by 2, contradicting optimality of S. This argument shows 
that i(tm) > m. If for t' > t^, < Tn, then we can obtain a more efficient 
solution by setting Si^t = for i > m and t <t' . 

Let S'i be the matrix consisting of the rows with index i > m. This 
matrix is a solution of the pebble game with n — m squares. Since this solution 
achieves the final configuration before tg, and by the choice of m, can use at 
most S — 1 pebbles. 

Since T{S) = T{Si) + T{S2) + T(S3), the desired inequality follows. 

To show that F(n, S) < mm(F(m, S) + F(m, S — 1) + F(n — m, S — 
1) I 1 < m < ra), we use the decomposition for the lower bound and show that 
it can be used constructively for an upper bound. Given Si, S2 and S3 with 
'^('^i) < S , S(S2) < 5* — 1 and SlSs) < S — 1, we construct a solution for the 
pebble game with n squares by composing the solutions Si. Let 4,8 and tg^i be 
defined for Si as we defined 4 and for S above (if m = 1 or ra — m = 1, 
set the corresponding = 4 + 1)- By shifting the solutions, we can assume 
t^ i + 2 = 4,3 and ie,2 = 4,3 + 2. Define 



'Sl;i,t 


for i < m 


and t < 4,1 + 1, 


*^3;«' — m,t 


for m < i 


< ra. 


'52;j,2te,2-t 


for i < m 


and 4,3 + 1 < 4 


1 


for i = m 


and 4,3 < t < te,3 + 1, 





otherwise 





Then 5 is a solution with F(S) = J2i satisfying the space bound. ■ 

The recursion of Theorem 2.1 can be used to compute tables of _F(ra, S). 
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The appendix contains some Mathematica expressions which compute F(n, S) 
dynamically. Tables 1 and 2 contain F(n, S) for small values of n and S . The 
proof of the theorem can be used to easily obtain the optimal strategy for each 
n and S . 

2.3 Properties of F{n,S) 

Some immediate observations are as follows: 

Observation 2.2 (i) F(n, S) = 2n - 1 for S > n. 

(ii) For S < S', F{n, S) > F{n, S'). 
(Hi) For n < n' , F{n, S) < F{n' , S). 

We next establish the maximum n for which F(n, S) < oo. 
Theorem 2.3 F(n,S) < oo iff n < 2^"^ 

Proof. The statement holds for S = 1. We proceed by induction on S and n. 
If ra < 2"^"^ and we let m = [ra/2j , then m < 2"^"^ and n — m< 2"^"^. It follows 
that F(m,S — 1), F(n — m, S — 1) and F(m,S) are finite by the induction 
hypothesis. The recursive expression for F(n, S) now implies that F(n, S) is 
finite. 

If ra > 2"^"^ then either m > 2"^"^ or ra — m > 2"^"^. Hence every term 
in the minimum of the recursion for F(n, S) is infinite. ■ 

To obtain explicit bounds on F(n, S) it is convenient to consider the 

function 



For ra > 1 and S > 1, let m(n,S) be the least minimizer of the recursive 
expression for F(n,S). Define m(0,S) = m(l,S) = 0. Note that m(n,S) < 
[ra/2j by the monotonicity properties of F(n, S). 

Lemma 2.4 . For m = m(n, S), ra > 1 and S > 2, 




if ra < 0, 

F(n + 1,5)- F(ra, S) if ra > and F(n + 1, 5) < oo. 



00 otherwise. 



A(ra-m,5-l) > A(m- 1,5) + A(m- 1,5- 1), 
A(ra-m- 1,5-1) < A(m, 5) + A(m, 5 - 1). 
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Table 1. Table of F(n, 5) for 1 < ra < 50 and 1 < 5 < 20. 
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CO 


CO 


CO 


CO 


CO 


CO 


481 


421 


381 


337 


291 


273 


271 


269 


267 


265 


263 


261 


259 


76 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


493 


429 


389 


345 


297 


277 


275 


273 


271 


269 


267 


265 


263 


77 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


505 


437 


397 


353 


305 


281 


279 


277 


275 


273 


271 


269 


267 


78 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


517 


445 


405 


361 


313 


285 


283 


281 


279 


277 


275 


273 


271 


79 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


529 


453 


413 


369 


321 


289 


287 


285 


283 


281 


279 


277 


275 


80 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


541 


461 


421 


377 


329 


295 


291 


289 


287 


285 


283 


281 


279 


81 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


553 


469 


429 


385 


337 


301 


295 


293 


291 


289 


287 


285 


283 


82 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


565 


477 


437 


393 


345 


307 


299 


297 


295 


293 


291 


289 


287 


83 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


579 


485 


445 


401 


353 


313 


303 


301 


299 


297 


295 


293 


291 


84 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


593 


493 


453 


409 


361 


319 


307 


305 


303 


301 


299 


297 


295 


85 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


607 


501 


461 


417 


369 


325 


311 


309 


307 


305 


303 


301 


299 


86 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


621 


509 


469 


425 


377 


331 


315 


313 


311 


309 


307 


305 


303 


87 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


635 


517 


477 


433 


385 


337 


319 


317 


315 


313 


311 


309 


307 


88 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


649 


525 


485 


441 


393 


343 


323 


321 


319 


317 


315 


313 


311 


89 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


663 


533 


493 


449 


401 


349 


327 


325 


323 


321 


319 


317 


315 


90 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


677 


541 


501 


457 


409 


357 


331 


329 


327 


325 


323 


321 


319 


91 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


691 


549 


509 


465 


417 


365 


335 


333 


331 


329 


327 


325 


323 


92 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


705 


557 


517 


473 


425 


373 


339 


337 


335 


333 


331 


329 


327 


93 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


721 


565 


525 


481 


433 


381 


345 


341 


339 


337 


335 


333 


331 


94 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


737 


575 


533 


489 


441 


389 


351 


345 


343 


341 


339 


337 


335 


95 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


753 


585 


541 


497 


449 


397 


357 


349 


347 


345 


343 


341 


339 


96 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


769 


595 


549 


505 


457 


405 


363 


353 


351 


349 


347 


345 


343 


97 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


785 


605 


557 


513 


465 


413 


369 


357 


355 


353 


351 


349 


347 


98 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


801 


615 


565 


521 


473 


421 


375 


361 


359 


357 


355 


353 


351 


99 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


817 


625 


573 


529 


481 


429 


381 


365 


363 


361 


359 


357 


355 


100 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


833 


635 


581 


537 


489 


437 


387 


369 


367 


365 


363 


361 


359 



Table 2. Table of F(n, S) for 51 < ra < 100 and 1 < 5 < 20. 
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Proof. Except for boundary cases, the result is obtained by comparing the 
expression in the minimum of the recursion of F(n,S) for m = m{n, S),m = 
■m(n, S) — 1 and m = m{n, S) -\- 1. The boundary cases are readily checked. ■ 

Lemma 2.5 The following hold: 

(i) A(ra, S) is nondecreasing in n. 

(ii) Let n > I and S > 2. If A(n - m(n, S), S - 1) < A(m(n,S),S) + 
A(m(ra, S), S —1), then m(n +1,5) = m(n, S). If not, then m(n +1,5) = 
m(ra, S) + 1. 

(iii) Let n > 1 and S > 2. Then either A(ra, S) = A(n — m(n, S), S — 1) or 
A(ra, S) = A(m(ra, S), S) + A(m(ra, S), S — 1), whichever is smaller. 

Proof. We prove the lemma by induction on S and n. For 5=1, A(ra, S) is 
nondecreasing in n by inspection. Let 5 > 1. We have m(0, 5) = m(l,S) = 
by definition. Note that m(2, S) = 1. A(0, 5) = and A(l, S) = 2 which gives 
the statement for ra = and n = 1. Let n > 2. By induction, A(n',S') is 
nondecreasing for S' = S — 1 and for S' = S and n' < n. Let m = m(n,S). 
Using Lemma 2.4 we get for m' < m 

A(ra + 1 - m'- 1,5- 1) > A(ra-m,5-l) 

> A(m- 1,5) + A(m- 1,5- 1) 

> A(m',5) + A(m',5- 1), 

which means that m' ^ m{n + 1,5). Similarly, for m' > m + 1, 

A(ra + 1 - m',5- 1) < A(ra-m- 1,5-1) 

< A(m,5) + A(m,5- 1) 

< A(m'- 1,5) + A(m'- 1,5- 1), 

so that m' ^ m{n + 1,5). This leaves only m{n +1,5) = m or m{n + 1,5) = 
m+1. The former holds iff A(n — m, 5—1) < A(m, 5) + A(m, 5—1). In this case, 
A(ra, 5) = A(ra — m, 5 — 1) and otherwise A(ra, 5) = A(m, 5) + A(m, 5 — 1). In 
the former case, by induction A(ra —1,5) = A(ra — 1 — m, 5 — 1) or A(ra —1,5) = 
A(m —1,5) + A(m — 1, 5 — 1), whichever is smaller, so that by monotonicity, 
A(ra,5) > A(ra-1,5). In the latter case, A(ra-1,5) = A(ra- 1 - (m- 1), 5 - 1) 
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or A(ra —1,5) = A(m — 2, 5) + A(m — 2, 5— 1), whichever is smaller, and again 
it is clear that A(ra, S) > A(n — 1, S). ■ 



Lemma 2.6 Let m = 171(11, S). Then 

A(n - m - 1,5 - 1) < A(n,S) < A(n-m,S -1), 
A{m-l,S) + A{m-l,S-l) < A{n,S) < A{m, S) + A{m, S - 1). 

Proof. This is an immediate consequence of Lemmas 2.4 and 2.5. ■ 

Lemma 2.7 A(n,S) is nonincreasing in S. 

Proof. The lemma follows from (i) of Lemma 2.5, Lemma 2.6 and by induc- 
tion. Let in = m(n, 5-1-1) and m' = m(n, S). We have 

A{n,S + l) = mm{A{n- m,S + l),A{m,S + 1) + A{m,S)) 
< mm(A(n — m, S), A(m, S) + A(m, S — 1)). 

If m' = m, then the last expression is A(ra, S) and we are done. If m' > m, then 
by Lemma 2.6, 

A{n,S) > A(m'- 1,5) + A(m'- 1,5- 1) 
> A(m,5) + A(m,5- 1). 

Similarly, for in' < in, 

A(n,5) > A(n-m'- 1,5- 1) 
> A(ra — in, 5). 

In either case, the desired inequality follows. ■ 

For A; > 0, let x^^s = min(ra | A(ra, 5) > 2^). This is well defined and 
nondecreasing in k. We have xq^s = 1, xi^s = S , and for A; > 0, Xk^i = 1. 

Lemma 2.8 For 5 > 2 and k > we have 

Xk+l,S > Xk+l,S-l + Xk,S-l, 
Xk+l,S < Xk+l,S-l + Xk,S, 

Xk+i,s < 2a;fc+i^5_i. 
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Proof. The proof is by induction on S and k. By Lemma 2.6, if either A(ra — 
m,S -1) < 2*^+1 or A(m, S - 1) < 2*^, then A(ra, S) < 2*^+^ This must occur 
if ra < Xk+i,s-i + Xk,s-i so that Xk+i,s > Xk+i,s-i + Xk,s-i- 

Suppose that x^^s > x^^i^s-i and n > 2xk^i^s-i- Since m = m{n, S) < 
[ra/2j , n — m > x^+i^s-i- If m = n — m = x^+i^s-i ■, then A(ra — 771,5— 1) > 
and A(m, 5) + A(m, 5 - 1) > A(m, 5 - 1) > 2*^+^ which gives A(ra, 5) > 2^+^. 
If m < a;fc_|_i^5_i , then n — m > x^+i^s-i and by Lemma 2.6, 2*^+^ < A(ra — m — 
1, 5 — 1) < A(ra, 5), so the last inequality follows. 

Suppose that x^^s < x^^i^s-i and n > Xk^s + 2;fc+i,5_i. If m = 
m(n,S) > Xk^s, then the second inequality of Lemma 2.6 and monotonicity 
of A imply A(ra, S) > 2*^"*"^. If m < Xk^s, then n — m > Xk^i^s-i and the first 
inequality of Lemma 2.6 can be applied to show that A(ra, S) > 2*^"*"^. Suppose 
that m = Xk^s and ra = x^+i^s-i- Then both A(ra — 777,5—1) > 2*^+^ and 
A(m, 5 — 1) > A(m, 5) > 2*^ and again A(ra, 5) > 2*^+^. This completes the 
proof of the lemma. ■ 



Lemma 2.9 For k > and S > 2, let 

k 

E 




mm 



Then x^l^g < Xk,s < a;^"5. 

Proof. It is straightforward to check that satisfies the recursion a;^'^-^ ^ = 
x^+i 5-1 + -^i's-i with the initial condition = 1. The function y^^s = 

i^~^k~^) satisfies the recursion yk+i,s = yk+i,S-i + yk,s with the correct initial 
condition. The bound 2"^"^ in the definition of a;^") comes from Theorem 2.3. 
The result follows by Lemma 2.8. ■ 

Lemma 2.9 cannot be used to obtain "nice" expressions for F(n, S) 
but yields lower and upper bounds by summing A(ra, S). 

Lemma 2.10 For 1 < k < S - 1 and S > 2 
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F(x^:ls) > 




Proof. The inequalities are obtained by summing A(n, S) and making use 
Lemma 2 

Y^fc (S+i-2\ 
l^i=0 \ i J 



of Lemma 2.9, ignoring the second case for a;^"^. Observe that ^ 



2.4 Asymptotic bounds 

For asymptotics where S = O(log(ra)), let 

f(l,S) = ^logF[[2^'\,s). 

All logarithms are base 2. The floor/ceiling notation will be omitted (or ab- 
sorbed into the 'o' and 'O' notation) for real numbers unless needed for clarity. 
The definition of f{-f,S) implies that F{n,log{n)/-f) = n^(7,losW/7)/7+o(i)_ 
For < 7 < 1, let 77(7) denote the usual entropy function 

H(j) = -7log(7) - (1 - 7)log(l - 7). 

Theorem 2.11 Let < 1/2. Then 

fiH{i),S) < 7 + i7(7) + o(l), 
/((l + 7)i?(7/(l + 7),5') > 7 + (l + 7)i?(7/(l + 7)) + o(l), 

where the o(l) are with respect to S 00. 

The proof of the theorem is below. First we make a few observations. 
For small 7, the two expressions in the bounds are the same to within 0(7^). 

H{-i) = -7log(7) + log(e)7 + 0(72) 

= 7(log(e/7)) + 0(7'), 
(1 + 7)1^(7/(1 + 7)) = -7log(7) + (l + 7)log(l + 7) 

= 7(log(e/7)) + 0(7')- 

So for small 7, the bounds are asymptotically tight. The upper bound is an 
asymptotic version of that given in [2] with a slightly better constant. More 
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specifically, let e = j^^^. Then 7 = 62"!/^ and H{^) = (1 + 7)1^(7/(1 + 7)) + 
0(7^) = + 0(7^). Substituting this in the theorem gives 

F (^n, ^2i/nog(n)) = „i+e+2-°^(i/^)+o„(i)^ 

where we used suffixes on the 'O' and 'o' notation to clarify the limit involved. 
Proof of Theorem 2.11. A consequence of the Stirling approximation of nl 
is that for a constant, a < 1/2, 

E-io (f) = o {O) = 2^(«)^(i+°(i)). 

Thus, for n = 2^W^, A(n,S) < 2"'^(i+°(i)). Using these observations and 
Lemma 2.10, one can estimate 

F{n,S) < (^^^^^2"'^^'+°^'^^ 

< 2(7+^^(7))5(l + o(l))_ 

Similarly, for a constant, 

x;"-o ("^^0 = o (('^^^5"'')) = 2(^+")^("/(^+"))'^(^+°(^)), 

so that for n = 2(i+"')^("'/(i+"')), A(n,S) > 2"'^(i+°(i)). This gives the second 
inequality. ■ 



2.5 Asymptotics of TS. 

The minimum time-space product is asymptotically determined by the next 
theorem. 

Theorem 2.12 

TSin) = 22^^(i+°(i))n. 
Proof. To show that TS(n) < 2Vi°sW(i+°(i)), let S = 2V^°s(^. Then 
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This implies that A(n, S) < 2^°sW(i+o(i)) and 

F(n,S) < 2^°s(^(i+°(i)), 



For the lower bound, we estimate F(n, S) for S = 2''v^^°^'^ from below. 
^ /5 + c^4^-l\ 

c^y\og{n),S - y cyiog(ra) J 

< (^5 + c^log(n)j 

< 2cdlog(n)(l + o(l)) ^ ^cd(l + o(l)) _ 

This implies that A(n, S) > 2^^^¥^(i+°(i)) and 

F(n,5)5 > 2^°s(^('^+d)(^+°(^))n(l - o(l)) 
> 22ViogW(i+o(i))^^ 

since rf + i > 2 for all > 0. We used the observation that the last term of the 
sum Yli=o asymptotically dominates it. ■ 



A Appendix 

A.l Mathematica expressions for dynamic evaluation of F{n, S) 

FBen[l,S_] := If [S >= 1, 1, Infinity]; 
FBen[n_,0] := Infinity; 

FBen[n_,S_] := FBen[n,S] = Module [{m>, 

Min [Table [FBen[m,S]+FBen[m,S-l]+FBen[n-m,S-l] ,{m, l,n-l>]] 

]; 

If the goal is to compute a table, then Lemma 2.5 can be used for 
efficient computation by keeping track of m(ra, S). 

FBenTable [nmax_ ,Smax_] := Module [{m,n,nt,S,T,Tt,tl,t2,t}, 
T = {Join[{l}, Table[Inf inity,{n,2,nmax}]]}; 
For[S=2, S<=Smax, S++, 
m=l ; 
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Tt = {1,3>; 

For[n=3, n<=nmax, n++, 

t l=Tt [ [m] ] +T [ [S- 1 ,m] ] +T [ [S-1 ,n-m] ] ; 
t2=Tt [ [m+1] ] +T [ [S-1 ,m+l] ] +T [ [S- 1 ,n-m- 1] ] ; 
If [t2<tl, 

m++; 

t = t2; 

J 

t = tl; 

]; 

Tt = Append [Tt,t] ; 
If [t==Inf inity , 

Tt = Join[Tt,Table[Inf inity, {nt,n+l,nmax}]] ; 

Break [] ; 
,Iull] 

]; 

T = Append [T,Tt] ; 

]; 

Return [T] ; 

]; 
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